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ABSTRACT 

This  study  is  a  continuation  of  the  investigations  in  the  use  of  automatic- 
throttle  compensation  systems.    The  previous  investigations  demonstrated 
the  capability  of  such  systems  in  alleviating  the  problems  of  aircraft  speed 
instability  in  the  carrier  landing  approach.     The  feasibility  of  adapting  the 
system  with  sensed-input  variables  of  angle-of -attack  and  normal  accelera- 
tion to  the  condition  of  field  landings  with  flare  was  investigated.    An  analy- 
sis  of  the  ground  effect  prevalent  in  fiefoiandings  was  performed  and  the 
significant  influences  on  aircraft  response  were  included  in  the  investiga- 
tion.    The  desired  airspeed  control  with  zero  rate  of  descent  at  touchdown 
was  obtained  by  incorporating  additional  commands  in  the  basic  auto-throttle 
device  during  the  flare  maneuver.     The  analog  computer  was  utilized  for 
problem  simulation. 

This  investigation  was  performed  byLCDR  Robert  L.   Elictv,   USN>,  at  the 
Uo  S.  Naval  Postgraduate  School,    Monterey,   California. 
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1.         Introduction. 

This  study  is  a  continuation  of  the  investigations  being  conducted  at 
the  U.  S.  Naval  Postgraduate  School  on  automatic-throttle  compensation 
systems.  The  previous  studies  conducted  by  Lieutenant  Commander  Go  R. 
Bell  [l]  ,  and  Lieutenant  R.  E.  Evans  and  Captain  R.  H.  Schuppe  [2]  were 
presented  in  1963  and  1964  respectively. 

As  shown  in  Ref.    1  and  Ref„   2  automatic-throttle  compensation  system 
utilizing  various  sensed-input  variables  are  capable  of  alleviating  the  speed 
instability  difficulties  of  modern  supersonic  aircraft  during  carrier  landing 
approaches.    One  such  system  utilizes  angle-of -attack  change  and  change 
in  acceleration  normal  to  glide  slope  for  commands  to  initiate  the  proper 
change  in  engine  thrust.    The  purpose  of  this  investigation  was  to  determine 
the  feasibility  of  adapting  this  system  to  the  condition  of  normal  field  land- 
ing which  is  terminated  with  a  flare  maneuver  to  the  touchdown.     It  is  pro- 
posed that  by  incorporating  additional  commands  to  the  basic  auto-throttle 
system  during  the  flare  maneuver,  a  decrease  in  airspeed  to  a  desired  min- 
imum at  touchdown  can  be  accomplished  simultaneously  with  zero  rate-of- 
descento 

The  linearized,  full  equations  of  motion  derived  from  small  perturba- 
tion theory  were  used  for  this  analysis.    The  coefficients  in  the  equations 
were  assumed  constant  except  for  those  directly  affeGted  by  the  aircraft 
flying  in  close  proximity  to  the  ground.  The  influence  of  ground  effect  on 
the  dynamic  stability  due  to  the  varying  coefficients  was  investigated  and 
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was  included  in  this  study. 

The  vehicle  used  in  this  and  the  previous  studies  is  the  Navy  swept 
wing  fighter,   F-8.    The  Dpnner  3100  Analog  Computer  with  associated  equip- 
ment was  utilized  to  simulate  the  aircraft  motion  and  maneuvering  response 
with  the  proposed  auto-throttle  device. 
2c         Basic  Automatic-throttle  (  cK  -  n)  System. 

The  type  of  automatic-throttle  controller  presently  used  in  the  Navy 
F-8  is  Specialties,   Inc.  Automatic-Power  Compensator,    (APC),  which  utilizes 
sensed-input  variables  of  angle-of-attack,  o<  ,  and  normal  acceleration,  n, 
for  engine  thrust  command.    A  block  diagram,  typical  of  this  system  is  shown 
in  Figure  1.     For  the  design  application  of  this  system,   constant  glide  slope 
and  airspeed  are  desired  throughout  the  approach  and  the  two  sensed  vari- 
ables of  instantenous  angle-of-attack  and  acceleration  normal  to  the  glide 
slope  in  "g"  units  are  the  engine  command  functions.    The  use  of  these 
variables  as  the  command  functions  for  airspeed  control  is  evident  from  the 
relationship: 

L    =    n'w    =    -^pV2S  CL 

n'W    =    j/?V2SCLo<o( 

n'     =      1   ^v2  S    Pt 

ill  =    Const.   (V2) 
From  this  equality,  and  increase  of  velocity  results  from  a  positive  normal 
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acceleration  increment  and  a  velocity  decrease  results  from  a  positive  angle- 
of -attack  increment.    Thus,  the  normal  acceleration  in  "g"  units  and  angle- 
of -attack  are  inputs  which  are  converted  into  the  thrust  command  by  the 
system  computer  utilizing  chosen  gain  and  time  factors.    The  resulting  change 
of  thrust  from  the  engine  acts  as  a  forcing  function  to  the  airframe  modes  of 
motion  appropriately  changing  the  aircraft  velocity.    The  APC  computer  in- 
directly controls  the  aircraft  velocity  through  the  equations: 
AT    =    K^A^    +    KnAn' 


since  AV   =    K-pAT 

therefore,        AV   =    K^Ao<  +    K^An' 

where  the  K  values  are  the  appropriate  shaping,  filtering  and  simulation 

networks  o 

The  transfer  function  for  the  thrust  command  as  specified  by  Specials 
ties,.  Inc.  is: 


AT      =    -, i 1 

c       (i  +  res)     (i+rQs) 


Kn 

n 


+— —  + ex 


1  +  7^  s    s        i  +  s  / 


where. 


i+rns 


7"     -  o.  2  sec  (throttle  servo  time  constant) 
s 

Tqi    -  0.5  sec  (angle  of  attack  filter  time  constant) 

7^     -  1.2  sec  (engine  time  lag) 

-7"     -  1.  0  sec  (normal  acceleration  filter  time  constant) 

K^-  570  lb/deg.   (angle-of-attack  displacement  gain) 
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K3    -    190  lii/deg-sec  (angle- of- attack  intergral  gain) 
K/    -    1,  100  lb/deg  (angle-of -attack  gain) 
Kn    -    19,  400  lb/"g"  (normal  acceleration  gain) 
As  evidenced  from  the  transfer  function,  the  effect  of  increasing  the 
instantaneous  angle-of-attack  results  in  an  opening  of  the  throttle  and  an 
increase  of  thrust,  and  a  decrease  in  thrust  is  obtained  for  a  decrease  in 
angle-of-attack.    An  increment  of  normal  acceleration  is  seen  as  a  negative 
element  in  the  transfer  function  with  the  effect  of  decreasing  the  thrust  for 
a  postive  "g"  perturbation  and  increasing  the  thrust  for  a  negative  "g"  per- 
turbation.   Ref.   1,   2,  and  3  support  the  validity  of  this  auto-throttle  device 
for  use  in  carrier  landing  approaches. 
3»         Discussion  of  Auto-throttle  System  in  Field  Landings. 

In  the  field  landing  phase,   it  is  desired  that  just  prior  to  touchdown 
the  vertical  kinetic  energy  be  dissipated  and  the  horizontal  kinetic  energy 
be  a  minimum  contingent  on  the  minimum  safe  airspeed.    Thus,   structural 
damage  and  wear  and  tear  on  aircraft  components  will  be  minimized  along 
with  the  advantage  of  using  shorter  runways.    From  the  pilots  point  of  view, 
these  results  are  attained  by  application  of  elevator  to  flare  the  aircraft  as 
necessary  to  decelerate  to  the  desired  touchdown  airspeed  with  zero  rate- 
of -descent. 

It  is  proposed  that  the  oc  -  n  auto-throttle  system  can  be  utilized  to 
accomplish  the  desired  airspeed  control  by  incorporating  a  supplementary 
circuit  in  the  designed  APC.    For  the  purpose  of  this  investigation,  it  was 
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assumed  that  the  F-8  aircraft  was  making  an  Instrument  Landing  System  (ILS) 
approach  with  a  four  degree  glide  slope  on  the  recommended  approach  speed 
of  1.  2  VgL  =139  Kt.  commencing  a  flare  at  100  ft.  above  the  ground,  approx- 
imately eight  seconds  prior  to  touchdown.    The  desired  airspeed  for  touch- 
down was  1.  05  VSL  =  123  Kt.  or  209  ft/sec. 

Since  it  is  mandatory  that  the  auto-throttle  system  be  unchanged  up 
to  the  flare  maneuver,  it  was  decided  that  the  proposed  supplementary  cir- 
cuit of  the  APC  be  separately  energized  at  this  point.    This  circuit  would 
supply  the  auto-throttle  device  with  additional  angle-of -attack  and/or  normal 
acceleration  signals  changing  the  error  signals  received.    Thus,  in  effect  a 
new  base  for  the  auto-throttle  system  is  created.    The  error  signal  changes 
shall  be  referred  to  as  "biasing"  throughout  the  remainder  of  this  study. 

The  method  of  applying  the  angle -of -attack  bias  was  determined  from 
the  equation  describing  the  change  in  pitching  moment  coefficient  from 
straight  to  rotational  flight: 

*Cm   =   Cm^*  +    Cmj)*n  +    Cmqq 

Assuming  that  the  change  in  moment  is  zero  and  neglecting  the  change  in 
moment  due  to  rate  of  pitch,   q,  the  resulting  equation  is: 

Cm 

^<*  «=  Const.    (A T/) 
Thus,   the  angle-of -attack  bias  was  experimentally  utilized  as  some  percent- 
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age  of  A  Tj . 

Since  it  has  been  previously  determined  that  the  average  pilot  pulls 
a  maximum  of  +  0.  2  g  above  the  1  g  steady  state  condition  during  the  flare 
maneuver    4],  the  normal  acceleration  bias  was  experimentally  utilized  as 
a  percentage  of  0.2  g. 
4o         Equations  of  Motion. 

The  basic  three  degree  of  freedom  longitudinal  equations  of  motion   . . 
based  on  body  axes  used  for  this  analysis  are  shown  below„    These  equations 
were  derived  from  linearized  theory  based  on  small  disturbances  from  a  ref- 
erence steady  state  flight  condition.    Full  derivation  and  definitions  are 
included  as  Appendix  I.     The  equations  are  shown  in  the  dimensional  form 
necessary  for  analog  simulation. 

X    FORCE: 

X„  Xa/  —  Xat 

T 


u*  -  %  u-^«  +gcos©e  =  -5jUt  A 


Z    FORCE: 


^U  -  Zw  Z  77       _ 

m  ml' 


mU0  m  niU0 


MOMENT  EQUATION: 


M„             MAc      >         Ma            ••        M„     .        Mat  Mt,^, 

JL  u  _  , — 2L  ex —  cx  +  0 2-  0  =  — =-i-AT  +  — 2-77 

R  R  R  RRR/ 


B  B  B  B  B  B 

The  values  of  the  dimensional  stability  derivatives  and  aircraft  parameters 
as  obtained  from  Ref.  3  are  included  as  Table  I  and  Table  II  respectively. 
5.         Ground  Effect. 

The  principal  effects  of  the  ground  on  a  nearby  aircraft  are  increases 
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of  the  slope  of  the  lift  curves  of  the  wing  and  tail  and  a  decrease  in  the 

downwash  at  the  tail.    These  changes  result  from  the  constraint  of  the 

vortex  patterns  at  and  behind  the  wing.     The  effects  may  be  approximated 

by  superimposing  an  image  vortex  system,  as  shown  below,   beneath  the 

ground  having  induced  velocities  of  necessary  magnitude  to  satisfy  the 

boundary  condition  of  no  flow  across  the  ground  plane   .    This  method  of 

approximation  utilizes  only  the  trailing  vortices  as  it  is  normally  considered 

that  the  bound  vortex  patterns  cancel  each  other.     The  trailing  image  vortex 

will  decrease  the  downwash  at  the  tail  and  decrease  the  downwash  in  the 

tail  and  wing  planes  which  increases  the  slopes  of  the  lift  curves. 

k Kb >| 


^V        ^ 


f^  REAL    _ 

y'  VORTEX 


G-ROUhJD 


IMAGE 

VORT  C  X 

For  a  two  dimensional  lifting  surface,  the  Kutta-Joukowski  theorem 
states  that  the  lift  is  equal  to  the  density  times  velocity  times  circulation 
or 

L   =yoV  F       lb/ft.  of  span 

^Andrews,   E.  J.,   Notes  on  Simplified  Ground  Effect  Analysis.    U.  S. 
Naval  Postgraduate  School,   Monterey,  California 
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Thus  the  local  lift  coefficient  is: 

ci  = 

/oV    P 

^/°\!2c 

mmm 

2T 

c   V 

For  an  airfoil  with  uniform  span  loading  with  the  associated  trailing  vortices 
a  distance  Kb  apart,   where  b  is  the  span,   the  induced  velocity  at  the  aircraft 
centerline  is: 

w    =    2 


r 


4  77  Kb/2 


Since  the  downwash,  €  ,  is  defined  as  -rr"i   then 


€    = 


r 


Mv  V  Kb/2 

C^  c 

4  7T  Kb/2 


Therefore,  the  change  in  downwash  with  respect  to  the  change  in  the  local 

lift  coefficient  is  proportional  to  — ^- —  .     Resolving  the  induced  velocities 

Kb/2 

from  the  image  vortex  system  into  the  vertical  plane  at  the  aircraft  centerline 
results  in: 


w 


€    = 


Kb/2. 


Z+(2d)2"U    f(%2    +  (2d)2U 


] 


C1   c.V  (*£-) 

4zrr(i^)2    +  (2d)2j 


d€      = 
dQ 


Kb/2 


4  7T 


(Kb),2    +  (2d)2 
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d€      ^,  Kb/2 

dC  <^)2    +  (2d)2 

Therefore,   in  proximity  of  the  ground  the  change  of  downwash  with  respect 
to  local  lift  coeificient  is  proportional  to  the  following: 


6C 


l Kb/2 


Kb  (Kb}2        ,    ,x2 

L/G      |_2  <2>      +(2d) 


By  setting  up  a  ratio  between  ^pr-  out  of  ground  effect  and  §f~  in  ground 
effect,  the  following  equality  is  obtained: 


( a  cJg  = 


i 


i  +  P)2 


4d 

Aerodynamic  texts  state  that  for  the  lifting  line  theory  of  finite  wings 
the  trailing  vortex  filaments  can  be  assumed  to  be  0.  8  of  the  span  apart. 
Substituting  this  value  for  K  in  the  above  equation  yields: 


(4*-) 

\  a  cJg 


fd€  \         1  +  0.4  (3p 


The  ratio  of  %%-  in  ground  effect  to  ^r  out  of  ground  effect  is  only  depen- 
dent  upon  wing  span  to  height  above  ground,   b/d.     Therefore,  a  single  curve 
utilizing  the  coordinates  of  the  ratio  of  downwash  versus  b/d  ratio  can  be 
constructed.     This  curve  is  shown  in  Figure  2. 

For  a  finite  aspect  ratio  wing,  the  lift  coefficient  is  given  as 
CL    =    2  7T  (  <*  -  0(;) 
where  o(;  is  the  induced  angle-of-attack  which  is  equal  to  — =-  for  ellip- 

*  7Tk 
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tical  lifting  surfaces.     Substituting  this  value  into  the  above  equation  and 
taking  the  derivative  with  respect  to   o<  gives: 


L  2   7T  A 


6  CK  A    +    2 

By  utilizing  the  same  procedure  asinthe  down  wash  effect,  the  following  ratio  of 
slope  of  the  lift  curves  im and  out  of  ground  effect  was  obtained: 


(a«7      A  +  2[tTT7T¥2J 


A   +    2 


d' 

Fora  particular  aircraft  with  a  wing  of  known  aspect  ratio,  a  plot  of  the  ratio 
of  the  slope  of  lift  curves  versus  b/d  can  be  constructed.    For  the  F-8  with 
an  aspect  ratio  of  3.  39,   such  a  plot  was  obtained  and  is  included  as  Figure 
3. 

Utilizing  the  curves  of  Figures  2  and  3,  the  coefficients  in  the  equa- 
tions of  motion  containing  stability  derivatives  dependent  on  downwash  or 
slope  of  lift  curves  of  wing  or  tail  were  varied.    The  stability  derivatives 
effected  were  X^  ,  Zw,   Mq(,  and  M^.     The  definitions  of  these  derivatives 
in  Appendix  I  were  rearranged  as  shown  in  Appendix  Iliinorder  to  use  the 
above  derivation.    Figures  4a,    4b,   4c,  and  4d  graphically  present  the  value 
of  these  derivatives  versus  altitude  above  ground.    The  values  of  the  coef- 
ficients containing  these  derivatives  were  then  obtained  and  are  presented 
in  Figures  5a,    5b,    5c,  and  5d. 

The  proximity  of  the  ground  on  an  aircraft  causes  a  lowering  of  the 
drag  polar,   reduced  drag  level  for  lift  attained.    This  effect  was  considered 
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offsetting  the  increased  drag  resulting  from  the  increase  of  angle-of-attack 
during  the  flare  maneuver.     Therefore,  the  drag  was  assumed  constant 
throughout  this  analysis. 
6.        Analog  Construction  for  Complete  Problem  Simulation., 

The  equations  of  motion  in  dimensional  form  rearranged  for  analog 
simulation  with  necessary  scaling  factor  notation  included  are: 
auii    =    -.0602auii   +    -^L.    ao(  <x     -    31.88aQ9 
+    .00145aATAT 

~  Zw  _ 

a^  <=x     =    -.00113anu   +    — —  cu,  o<     +   a.0 

u  m  9 

-    .0595 


^V 


a..0      =    1.89x10  4aSlu   +    — —  a^,  «x     +    a^** 

0  u  B      ^  B 

-    .  3396a.  0     -    .455xl0"5a    ^AT  -  2.  504a„7? 
q  AT  / 

The  scaling  factors  are  included  in  Table  III  and  the  analognraechahiEfiition 
of  the  equations  is  shown  in  Figure  6. 

The  automatic -throttle  transfer  function  given  by  Specialties,  Inc. 
was  converted  into  analog  format  as  shown  in  Appendix  HI  with  the  scaling 
factors  included  in  Table  III  and  then  mechanized.     The  analog  mechaniza- 
tion is  shown  in  Figure  7. 

For  the  complete  problem  simulation,  auxiliary  equations  for  generat- 
ing actual  altitude  during  the  approach,  glide  slope  and  acceleration  normal 
to  the  glide  slope  in  "g"  units  were  derived  as  shown  in  Appendix  IV  with 
scaling  factors  also    ncluded  in  Table  III.      The  equations 
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were  mechanized  as  represented  in  Figure  8.     These  equations  in  analog 
format  are: 

Altitude: 

aHB   =  aHoH0    -  X0(md)J*uJj9  +    KiJ^+^u)J*t* 
Glide  slope: 

^^   =   ae9    -   a^ 
Acceleration: 

a'n  "f  =  "g  (au0,XJD+  au^)^e  S   -jfe  «  ) 

To  perform  the  flare  maneuver,  an  analog  mechanization  for  elevator 
control  was  constructed  by  applying  a  steady  voltage  through  a  potentiometer 
to  an  integrator  amplifier  circuit.     Thus,   by  adjusting  the  potentiomenter, 
application  of  step  input  or  variable  slope  ramp  input  of  elevator  to  the  air- 
frame could  be  obtained.    The  elevator  circuit  was  further  controlled  by  an 
off/on  function  switch  built  into  the  Donner  computer. 
7.         Experimental  Methods  and  Results. 

7.  1     Ground  Effect  Investigation. 

Due  to  an  insufficient  number  of  available  amplifiers  on  the 
Donner  Analog  Computer  for  the  simulation  required  in  this  study,   it  was  im- 
possible to  include  all  the  ground  effect  variables.    Therefore,   it  was  necess- 
ary to  evaluate  the  influence  of  each  variable  on  the  dynamic  stability  of 
the  basic  airframe  inorder  that  a  simplification  could  be  obtained.     Prior  to 
experimentation,   it  was  concluded  that  ground  effect  acting  on  the  derivative 
M  ^  would  be  the  most  influential  and  therefore  would  be  retained. 
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A  dynamic  stability  analysis  on  the  airframe  equations  of  motion  with- 
out  forcing  functions  of  thrust  and  elevator  was  performed.     Using  initial 
condition  of  u(0)  =  25  ft.  /sec.  for  the  phugoid  determination  and  an  initial 
condition  of  o<  (Q)  =  5  degrees  for  the  short  period  mode  of  motion  determina- 
tion,  the  following  analog  runs  were  obtained: 

(1)  Time  histories  of  ot  and  u  perturbations  out  of  ground  effect  in- 
cluded as  Figure  9. 

(2)  Time  histories  of  ok  and  u  perturbations  with  all  ground  effect 
variables  included  for  heights  above  ground  corresponding  to  b/d  ratios 
1.0/   2.0/   3.0,  and  4.  0.     Results  are  included  as  Figure  10. 

(3)  Time  histories  of  ex  and  u  perturbations  neglecting  ground  effect 
on  derivatives  X^  and  Zw  for  the  same  altitudes  as  in  (2).  The  analog  run 
outs  are  included  in  Figure  1 1 . 

(4)  Time  histories  of  o<  and  u  perturbations  neglecting  ground  effect 
on  derivatives  X^  ,  Zw,  and  M£  for  the  same  altitudes  as  before.    The 
analog  representations  are  included  in  Figure  12. 

Analysis  of  Figure  9  determined  that  the  phugoid  period  is  34  seconds 
and  that  of  the  short  period  mode  is  six  seconds.     Times  to  damp  to  h  ampli- 
tude are  approximately  30  seconds  and  2  seconds  respectively.    With  the 
condition  of  all  variables  included,    the  maximum  deviation  from  the  no 
ground  effect  analysis  in  the  phugoid  and  short  period  modes  was   0.  3  sec- 
onds;-       There  was  no  noticeable  change  in  the  times  to  damp    *    ampli- 
tude.     Analysis  of  the  analog  runs  when  holding  X^    and  Zw  constant 
showed  no  noticeable  change  of  period  or  time  to  damp  to  -  amplitude  of 
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either  mode  of  motion  from  those  obtained  with  all  derivatives  variable. 
When  the  derivatives  M^,    Zw,  and  X^  were  held  constant,  the  phugoid 
mode  period  was  0.  2  seconds  longer  and  the  short  period's  0.  1  seconds 
shorter  then  obtained  when  varying  all  the  derivatives.    The  times  to  damp 
to  y  amplitude  were  basically  unchanged. 

Since  the  relative  influence  of  the  variable  derivatives  X^ ,  Z^,  and 
M&   on  the  dynamic  stability  was  inconsequential,  the  decision  was  made 
to  conduct  the  investigation  of  this  study  with  these  three  derivatives  in- 
variant to  ground  effect.    The  derivative  M©<   which  was  retained  as  a 
variable  dependent  on  altitude  was  constructed  on  a  function  generator  for 
analog  use.    The  analog  mechanization  is  shown  in  Figure  13. 

7.2     Design  Auto-throttle  System  Operation. 

The  normal  auto-throttle  system  with  a^nd  without  ground  effect 
was  investigated  to  insure  proper  analog  mechanization  prior  to  commencing 
experimentation  for  the  proposed  system  during  the  flare  maneuver.    Analog 
runs  with  auto-throttle  operating  for  initial  conditions  of  horizontal  gusts, 
u(0)  =  +5  knots,  vertical  gusts,  ex  (0)  =  5°,  and  pitch  angle  9(0)  =  +  4°. 

Figure  14a  shows  the  response  obtained  without  ground  effect  for  an 
initial  horizontal  gust.    The  velocity  perturbation  is  oscillatory  damped  with 
an  initial  overshoot  approximately  45%  of  u(0)  and  was  without  a  second  over- 
shoot.    Figure  14b  is  the  response  for  an  initial  horizontal  gust  with  ground 
effect  and  agrees  closely  with  that  of  no  ground  effect.     The  thrust  applica- 
tion is  oscillatory  in  both  cases  with  a  damped  wave  motion  with  a  peak4  Is 
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1200  lb. 

Figures  15a  and  15b  present  the  response  obtained  for  an  initial  verti- 
cal gust  with  and  without  ground  effect  respectively„     The  oscillatory  veloc- 
ity perturbation  is  rapidly  damped  with  an  oscillatory  throttle  application 
with  a  peak  4  T  =  6,  300  lb. 

The  analog  results  for  the  initial  pitch  angle  condition  without  ground 
effect  is  presented  in  Figure  16.     The  velocity  perturbation  is  postively 
damped  with  a  slight  oscillatory  tendency  which  is  always  negative.     The 
thrust  application  is  smoother  than  in  the  other  initial  conditions  and  has  a 
peak 4  T  =  600  lb. 

The  results  obtained  for  the  design  application  investigated  correspond 
to  the  results  obtained  in  Ref.    1  and  are  as  expected.     The  analog  mechaniz- 
ation was  therefore  considered  correct.     The  numerical  forms  of  all  the  equa- 
tions are  as  in  Appendices  I  -  IV  and  the  analog  potentiometer  settings  is 
included  as  Appendix  V0 

7  =  3     Elevator  Angle  for  Flare  Maneuver. 

The  elevator  deflection  required  to  perform  the  flare  maneuver  was 
approximated  from  the  equation  given  in  Ref.   5  relating  the  elevator  angle 
per  "g"  for  a  constant  acceleration  maneuver.     This  equation  is: 

CT    + 


*V 


'mo<        L    T      \j    ^mq    ^L^ 


n  "  l  CL*  Cm7?     "    Cm^  CWj 

Substituting  in  the  required  values  for  the  non-dimensional  derivatives  and 

an  average  normal  acceleration  of  1.  1  g,  the  resulting  elevator  deflection 
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was  approximately  -2.7  degrees. 

Since  a  pilot  will  normally  pull  a  maximum  of  +  0.  2  g  above  the  steady 
state  1  g  during  flare,  the  rate  of  elevator  application  was  obtained  by  ad- 
justing the  potentiometer  in  the  elevator  circuit  until  a  peak  "g"  perturbation 
of  0o  2  g  was  obtained. 

Thus,   elevator  was  applied  at  the  required  rate  until  a  total  of  -  2.7° 
was  obtained  or  touchdown  occurred. 

7o  4     Determination  of  Auto-Throttle  Supplementary  Circuit. 

In  order  to  obtain  an  insight  as  to  what  commands  the  supple- 
mentary circuit  shall  supply  the  design  auto-throttle  system,  analog  runs  of 
a  flare  maneuver  with  constant  thrust  and  with  the  design  system  operating 
were  performed.     The  constant  thrust  response  is  shown  in  Figure  17.     The 
velocity  is  seen  to  decrease  18  ft. /sec.  during  the  flare  maneuver.    The 
response  obtained  while  the  auto-throttle  was  operating  is  shown  in  Figure 
18.    The  velocity  decreases  initially  then  attempts  to  return  to  the  beginning 
steady  state  velocity  as  a  result  of  an  increase  of  thrust. 

The  above  results  indicated  that  the  additional  commands  prescribed 
for  the  auto-throttle  supplementary  circuit  must  effectively  cancel  the  com- 
mands that  caused  the  thrust  to  increase. 

The  analog  mechanization  of  the  two  biasing  commands,  as  discussed 

in  Section  3,   is  shown  in  Figure  19.     To  simulate  real  system  operation, 

the  same  gains  were  used  in  the  angle -of -attack,  bias,   CX   ,  and  normal  accel- 

B 

eration  bias,   gB,  as  in  the  original  APC  system;  however,  an  additional 
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potentiometer  was  inserted  in  each  circuit  inorder  to  control  the  amount  of 
©<B  or  gB  sent  to  the  auto-throttle  block.     The  biasing  command  could  then 
be  comprised  of  both  c*g  and  gB  or  either  one  of  them. 

Analog  experimentation  for  determining  the  bias  commands  was  per- 
formed as  follows: 

(1)  Biasing  only  the  normal  acceleration  loop  through  a  range  of 
0.2  g  to  0.  01  g. 

(2)  Biasing  only  the  c*  loop  through  the  range  of  0.  01  AT)  to  1.0  LT). 

(3)  Biasing  both  the  c*  and  "g"  loops  through  their  respective  ranges. 

The  responses  for  the  normal  acceleration  biasing  were  initial  decreas- 
es in      airspeed  to  a  minimum  value  followed  by  an  increase.     The  minimum 
velocity  was  never  less  than  224  ft. /sec.  regardless  of  the  amount  of  gB. 

as  an  increase  of  thrust  was  always  obtained.     The  aircraft's  motion  was 
a  pull  up  with  a  variable  minimum  altitude  dependent  on  amount  of  gB.    A 
representative  record  of  the  typical  response  obtained  is  included  as  Figure 
20  which  is  the  run  with  gB  =  -  0.  02  g. 

The  typical  response  obtained  when  biasing  only  the  angle-of -attack 
loop  was  a  smooth  decrease  in  velocity  to  a  minimum  followed  by  a  slight 
increase  to  a  new  steady  state  value.     The  angle -of -attack  increased  rapidly 
to  a  peak  angle  then  decreased  to  a  damped  perturbation  about  some  median 
angle.     The  thrust  was  decreased  initially  by  a  maximum  of  three  to  to 

four  thousand  pounds  depending  on  the  amount  of  bias  and  then  began  increas' 
ing  as  the     flare  continued.    As  o<B  was  increased,  a  greater  decrease  in 
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velocity  and  higher  peak,  angle-of -attack  were  attained.    With  the  described 
elevator  application,  an    o<g  of  0.  80277  resulted  in  a  smooth  velocity  de- 
crease of  25  ft. /sec.  to  the  desired  209  ft. /sec.  at  the  touchdown  point. 
Touchdown  occurred  seven  seconds  after  flare  was  commenced  but  the  ver- 
tical kinetic  energy  was  not  completely  dissipated  as  a  slight  vertical  veloc- 
ity component  was  evidenced  in  the  altitude  read  out.     The  peak  angle-of  - 
attack  =  3.  8°  occurred  at  an  altitude  of  20  ft.  and  the  maximum  thrust  de- 
crease was  3000  lb.     This  presentation  is  included  as  Figure  21. 

The  results  obtained  during  the  analog  runs  while  biasing  both  the 
angle-of -attack  and  normal  acceleration  loops  indicated  that  the  airspeed 
response  could  be  more  advantageously  controlled.    Also,  the  magnitude 
changes  in  thrust  were  less  than  those  obtained  in  the  c<g  condition.    The 
peak  angles -of -attack  were  also  decreased  along  with  their  rate  of  increase. 
The  airspeed  increased  slightly  as  the  flare  was  commenced  then  smoothly 
decreased  to  a  minimum  value.     This  minimum  velocity  was  followed  by  a 
damped  perturbation  about  a  velocity  approximately  four  ft. /sec.  greater 
than  the  minimum.    With  the  elevator  application  as  previously  stated,  an 
optimum  condition  was  obtained  with  o<g  =  .  80577  and  gB  =-.  15g.    This 
bias  condition  resulted  in  the  smoothest  decrease  in  velocity  to  a  minimum 
209  ft. /sec.  that  occurred  exactly  at  the  touchdown  point.    The  flare  was  a 
smooth  manuever  with  essentially  all  the  vertical  kinetic  energy  being  dis- 
sipated.    The  peak  angle-of -attack  =  3.  3°  occurred  6.  5  sec.  after  flare 
commenced  and  3.  2  sec.   prior  to  touchdown.     The  thrust  increased  initially 
to  a  peak  of  +  400  lb.  then  decreased  to  a  minimum  -  2200  lb.  and  finally 


increased  to  a  positive  1000  lb.  above  the  beginning  steady  state  value. 
8.         Summary  of  Analog  Results  and  Recommendations. 

The  analog  experimentation  disclosed  that  a  desired  reduction  of 
velocity  could  be  obtained  during  the  flare  maneuver.    Two  different  supple- 
mentary circuits  to  the  design  auto-throttle  produced  a  required  touchdown 
velocity  with  a  simplified  elevator  application.     These  circuits  were  com- 
prised of  (a)  biasing  only  the  angle-of -attack  loop  by  a  percentage  of  the 
ramp  elevator  application,  and  (b)  biasing  both  the  angle -of -attack  and 
normal  acceleration  loops  with  a  percentage  of  the  elevator  and  a  constant 
value  of  "g"  respectively. 

The  o<b  produced  a  smoother  throttle  response;  however,   it  had  a 
greater  magnitude  change  of  thrust.    This  bias  also  caused  a  more  rapid 
o(  change  and  a  greater  peak  o<  while  not  dissipating  all  the  vertical  kinetic 
energy.     The  energy  could  be  dissipated  by  changing  the  elevator  deflection 
but  the  desired  airspeed  response  could  not  be  attained.    The  circuit  com- 
prised of  biasing  both  loops  had  a  less  smooth  thrust  response  with  a  very 
slight  velocity  increase  initially.     The  flare  maneuver  was  smoother  with 
essentially  ail  the  vertical  kinetic  energy  dissipated.     The  rate  of  c*  in- 
crease was  improved  and  a  less  peak  <X  was  obtained. 

In  general,  the  supplementary  circuit  entailing  biasing  both  auto- 
throttle  input  loops  produces  better  aircraft  response  even  though  the  thrust 
application  was  not  as  smooth.     Further  investigation  to  obtain  an  optimum 
supplementary  circuit  is  required  when  all  the  ground  effect  variables  are 
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included  along  with  the  changes  in  the  lift  and  drag  coefficients  resulting 
from  the  increased  angle-of -attack  during  the  flare..   An  analysis  in  system 
operation  with  gusts  during  the  flare  maneuver  and  the  associated  elevator 
application  is  also  required. 

It  is  recommended  that  the  required  further  analyses  be  performed  on 
a  hybrid  simulator  consisting  of  analog  and  digital  computers.    These  addi- 
tional analyses  could  not  be  performed  in  this  study  as  the  Donner  3100 
Analog  Computer  was  completely  saturated. 
9.         Conclusions  and  Acknowledgements. 

1.  The  use  of  the  auto-throttle  device  with  sensed  inputs  of  angle- 
of -attack  and  acceleration  normal  to  the  glide  slope  for  field  landings  with 
a  flare  maneuver  is  determined  feasible. 

2.  The  throttle  could  be  controlled  to  produce  a  reduction  of  airspeed 
during  the  flare  to  a  desired  minimum  at  touchdown  together  with  zero  rate 

of  descent, 

3.  The  desired  results  were  accomplished  by  incorporating  a  sup- 
plementary circuit,   energized  during  the  flare  maneuver,   in  the  design  auto- 
throttle  system.    This  circuit  would  be  comprised  of  an  angle-of -attack  rate 
command  to  a  maximum  and  a  constant  signal  command  of  normal  acceleration. 

4.  As  might  be  expected,    system  functioning  appeared  to  be  critical- 
ly dependent  upon  the  rate  of  application  of  elevator  during  the  flare. 

The  writer  wishes  to  express  his  gratitude  to  Professor  E.  J.  Andrews, 
Department  of  Aeronautics,   U.  S.  Naval  Postgraduate  School,   for  his  tech- 
nical assistance  and  encouragement  in  the  preparation  of  this  study. 
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TABLE    I 
LONGITUDINAL  STABILITY  DERIVATIVE  VALUES 


*u 

"41-14        ft. /sec. 

Xw 

= 

n    71                   lb. 

y  .  /  1           "           '   ;    ' 

ft. /sec. 

X<* 

= 

-2272.14       lb. 

X^T 

= 

n   oo               lb. 

u.  yy         ,      , 

ft. /sec. 

Zu 

= 

i  r>  i       /-n                        ID. 

ft. /sec. 

zw 

= 

oni      pq                      lb. 

ft„/sec. 

z7 

= 

-  9522.04       lb. 

Mu 

= 

18.15       lb. -sec. 

M<5< 

= 

-  3940.41       ft. -lb. -sec. 

Met 

= 

-  109,  377.45       ft. -lb. 

Mq 

= 

-  32,  602.94      ft. -lb. -sec 

M^r 

= 

-  0.437       lb. -sec. 

M^ 

= 

240,342,91     ft. -lb. 
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TABLE    II 
F-8  AIRCRAFT  PARAMETERS 


UQ   =    234  ft. /sec.  CD      =0.0 


e.g.     =    24%  MAC  CL        =    2.598 


W   =    22,000  1b.  CL        =    0.3897 

CX0    =    8.1°  C_      =    0.738   X    10"2 

o  mu 


0n    =    4.1°  Cm        ~    -  0.38 

o  m^ 


Sw    =    375  ftz.  Cm.      =    -  0.55 


ST    =    93.4  ft2.  CLo,T    =    1.622 


b   =    35.666  ft2.  h    =    0. 24 


lt    =    14.081ft.  hnwb   =    0.293 


c    =    11.78  ft.  Ct         l    =    2.387 

CL    =    0.87  -4^-    =    0.4772 

QO< 

CD    =    0.197  A    =    3.39 


Cn        =    0.963 
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TABLE    III 
ANALOG  SCALING  FACTORS 
Maximum  Voltage        100  volts 

u:         100    v.  =  50  ft.  /sec.  ;     au  =  0.  5  ft./sec. -v. 
u:         100  v.  =  50  ft./sec2.     ;     au  =  0.  5  ft./sec2.-v„ 
oc:       100  v.  =  10°  =  .  1745  rad.     ;     a^    =  0.  l°/v.  =  .  001745  rad/v. 
<*:       100  v.  =10°/sec.  -  .  1745  S&     ;     a<S<  =  .  1  ||a_ 
0:         100  v.  =  20°  =  .  349  rad.     ;     aQ  =  .  2°/v.  =  .  00349  rad/v. 
9:         100  v.  =  20°/sec.  =  .  349  ^r    ;     a;  =  .  2°/sec-v.  =  .  00349  jjft&  ■ 

act.        O  S6C~V« 

5:    100  v.  =  20°/sec2  =  .  349  ^2  ;  ax-  =  .  2°/sec2-v.  -  .  00349  &&-2 

sec     b  sec  -v. 

ff  :   100  v.  =  20°  =  .  349  rad.  ;  a  ^  =  .  2°/v.  =  .  00349  rad/v. 

n:    100  v.  =  2  g  ;  an  =  .  02  g/v. 

UQ:   100  v.  =  300  ft./sec.  ;  ay  =  3.  0  ft./sec. -v. 

o 

AT:   100  v.  =  8000  lb.  ;  a^T  =  80  lb. /v. 

77:   100  v.  =  28.65°  =  .5  rad.  ;  a^  =  28.  65°/v.  =  .  005  rad./v. 

^  100  v.  =  2.  0  /sec2.  ;  ^2L  =  t  02  /sec2-v. 

H:    100  v.  =  300  ft.  ;  aH  =  3.  0  ft. /v. 
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FIGURE  1 
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FIGURE  2 


GROUND  EFFECT  CSZ 
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FIGURE  3 
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FIGURE  4a 
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FIGURE  4b 
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FIGURE  4c 
STABILITY  DERIVATIVE  M  & 
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FIGURE  4d 


STABILITY  DERIVATIVE  Mc* 
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FIGURE  5a 
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FIGURE  5b 
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FIGURE  5d 
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FIGURE  6 


ANALOG  MECHANIZATION"  OF 
EQUATIONS  OF  MOTION 
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FIGURE  8 

ANALOG  MECHANIZATION"  OF 
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FIGURE   9 
AIRFRAME  STA3ILITY  ANALYSIS 
110  GROUND  EFFECT 
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FIGURE  10 
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FIGURE  11 

airframe;  stability  analysis 

DERIVATIVES  Xo<  ,   Moc   VARIABLE 


PPIUGOID  MODE 


SHORT  PERIOD  MODE 
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FIGURE  12 
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FIGURE  13 

ANALOG  MECHANIZATION  OF 
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FIGURE  14a 
AUTO-THROTTLE  FOR  u(O) =  5  kt. 

WITHOUT  GROUND  EFFECT 
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FIGURE  14b 
AUTO-THROTTLE  FOR  u(O)^  5  kt. 

WITH  GROUND  EFFECT 
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FIGURE  15b 

AUTO-TKROTTLB  FOR  <*(o)  -  ,5° 

WITH  GROUND  EFFECT 
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FIGURE   16 
AUTO-THROTTLE  FOR  0  (0)  *  4° 
WITHOUT  GROUND  EFFECT 
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FIGURE   17 


RESPONSE  DURING  FLARE 
WITH  CONSTANT  THRUST 
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APPENDIX    I 
EQUATIONS  OF  MOTION 
The  basic  three  degrees  of  freedom  aircraft  longitudinal  equations 
of  motion  as  derived  from  small  perturbation  theory. based  on  body  axes  are: 
X  Force: 

[mfj-   Xu]u    -    wXw   +  [mgcosG    -   Xqfje    =   X(t) 
Z  Force: 

-  uZu    +   |"m|-  -    Zw|w    -    f(mU0    +    ZqJ^Je    +    mg  sin©=  Z(t) 
Moment  Equation: 

-  uMu    -    wMw    -    wMw    «i  0B    -    9M-q    =    M(t) 
Rearranging,   substituting  in  forcing  functions  of  thrust  and  elevator,  and 
neglecting  minor  terms  Xq,  X^  ,   Z  ,   sin©,  and  Z^    the  equations  become: 


X  Force: 


U    -    -#  U    -    —  «    +    g  cos®  0    =    ~t~   A  T 


where  o<   = 


m  ~ 
w 


Z  Force: 


mUo      u  m  ,  mU0     / 


Z„u    +    OC    - 


Moment  Equation: 


B  B  B  B  B  B  * 
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Figure  defining  signs: 


u. 


Definitions  of  Stability  Derivatives: 

\i    =    7°uosw(cD   +    cdJ 

Xw  ■  i/>u0sw(cL  -  Co) 


X<X  =    U^ 


w 


\dT 


=    cos 


<«.+  £> 


Zu    =    -/°UoSw(CL    +    C,) 


u 


Zw    -    -f/oUoSwCC^    +    CD) 

*7    =    -^UoSwCL^ 

Mu    =  /<>UoSw^Cm    +    iUoCmu) 


M 


o< 


M, 


"    Uo("  2ST1TCLo<TarSr) 

=    U0MW    =    Up^/oU^  Cmo( ) 


Mq    =   -I  V'UqSwC    Cm  ) 


Kj  I   5    -  0.  437 


M 


7    =   ^UoS  *   Cm^  ) 
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Substituting  the  values  given  in  Table  I: 
X  Force: 

u   +    .0602u   +    3.325  c*  +    31.88  0    =    .001454T 
Z  Force: 

.00113  u    +    ck    +    .4271  ck     -    0    =    -.059577 
Moment  Equation: 

-  1.  89    x    10"4  u   +    .  041  cx  +    1.  139  o<  +    0    +    .  3396  0 
=    -  .455  x  10"54T    -    2.504  77 
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APPENDIX    II 


STABILITY  DERIVATIVES  FOR  GROUND  EFFECT 


Xo<-* 


Z  /°  Uo    SW  CL     ( 


/    - 


2    C^    ) 
&7T  A    I 


~k  /°Uo    S^C*. 


I     - 


&7T 

Z  „         (C^U 


e  tta  *"«   (C^  ) 


*w' 


w 


'  -7/°4,Sw  (cLm   *  cB) 


s  "  2  /°  Uo  $»/ 


'« frj      c 


U*M„  '  "c 


i* 


/oUcS^c   C^ 


] 


Cm     =  CL  (h  -  hn        )  -   VH  Cu        (l  -    4^~  ) 


=    C, 


(Ct^b)o 


w 


-  CT -{h  ■ h- ) 


-  Hi  CL. 


a«  Iftd-uHsJ*  (z&l 


Ha: 


CX 


M*    = 


=  U0M 


-J-         r      c*  r  $-€ 


ex 
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M. 


"2  /°°' 


=  -3740.+/ 


(c^)g    (^cjg.    (d-c^G 
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APPENDIX    III 


ANALOG  FORM  OF  AUTO-THROTTLE  EQUATIONS 


AT    = 


(i+7-gS)  (i+rss) 


Kc* 


1 


(1  +  1.2S)    (1  +  ,2S) 


1  +  ^S 


570 
1  +  .5S 


S 


Kf 


1  +  S 


o< 


190      1100   i~ 


1  +7"  S  • 

n 


19,400 
~~~ — —  n 
1  +  S 


Equation  without  time  elements 

AT    =    570(57. 3)©<   +    190(57.  Z)fcK   +    (1100)(57.  3)°<    -    19, 400n 
a^T^f  =    32,661a0<o<     +    10, 887  ^/o?  +    63,000^3!  -    19,  400  ann 

AT    =    .71  «    +    .237/o<    +     1.371  «     -    4.85  IT 
where: 

Ko<    =    570  1b./deg.     =    32,  661  lb./rad. 
Kj   =    190  lb./deg-sec.    =    10,  887  lb./rad-sec. 
Ky     =    11,  000  lb./deg.     =    63,  000  lb./rad. 
Kn     =    19,400  lb./g 
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APPENDIX    IV 


ANALOG  DERIVATION  OF  AUXILIARY  EQUATIONS 


Acceleration: 


n    =    h  (normal  to  fit.   path)    =    V 


Total 


=    (UQ    +    u)    ^-(0    -  c*  ) 

=  (u0  +  u)  (e  -  <*  ). 

a„n   =    (anU0   +    auu)    0^9        a^  «  ) 


'n1 


9 


.  02  n    =    (3.  0    U0    +    .  5  u)    (.  00349  9    -    .  001745  ©0 


^-("g"  units)    =   ^qT$o   +    .  167  u)    (.524    9    -    .  262  «  ) 


Glide  Slope: 

From  Appendix  I; 


8  =    9    -  <=* 


a*  #  = 


ae9    -    a«  ex 


1<X 


.00349  5    =    .00349  9    -    . 001745  ex 


er  =   9    -    .5  c* 


Total  Altitude: 


HQ    -    Initial  altitude    =    300  ft. 


H, 
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H    =    HQ    -    UQ(t)  sin  tf0    +    h 


=    Ho    - 

^(rad)/uo   +J/h 

=    Ho    - 

%/u0  +J/vT  i    . 

-    Ho    - 

.07_/uo  +  (u0  +  u)/v 

aHH    =    aH0Ho    "    -07/aUouo   +    <aU0Uo   +    auu)/a^ 
3. OH    =    3.0HQ    -    .07^3.0^   +    (3.0UQ   +    .  Su)  J^ .  00349  ^ 

H  =    i      -    .07/u_    +    (U.L  /-fo).  167  u)(.  0087)/-^- 
°  J        °  °  V     2.5 
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APPENDIX    V 


ANALOG  POTENTIOMETER  SETTINGS 


Potentiometer 

Coefficient  of  (Equat.) 

Setting 

1 

u  CuJ 

.0602 

2 

e  (u) 

.2252 

3 

<x(u) 

.0116 

4 

oC(<^) 

o4271 

5 

u  (<*) 

.3238 

6 

u  (e) 

o0271 

7 

«  (e) 

,0205 

8 

0  (e) 

.3396 

9 

o<  (tf) 

.5 

10 

VT(h) 

.87 

11 

tf(h) 

c5 

12 

o<  (KL    - 

-    T) 

.1371 

13 

n(Kn    - 

T) 

.485 

14 

o<  (Ky 

-    T) 

.237 

15 

©((K^ 

-    T) 

.71 

16 

SB 

variable 

17 

Y0(h) 

.093 

21 

77(e) 

.358 

22 

77  (<x) 

,1705 

23 

u0  (vT) 

.78 
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Potentiometer 

Coefficient  of  (Equat. ) 

Setting 

24 

u  (VT) 

.167 

25 

u0  (*/u0) 

.0275 

26 

H0  (H) 

1.0 

28 

0  (A  n) 

o524 

29 

c<  (4n) 

o262 

30 

g  (An) 

.311 

31 

<*B(K1    " 

T) 

.1371 

32 

^B^d     ~ 

T) 

.237 

33 

°<b(kcx     " 

T) 

.71 

34 

Funct.  Gen, 

.1 

35 

AT(u) 

.232 

36 

a  t  (e) 

.1043 

37 

°<B 

variable 

38 

T) 

variable 
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